Habitat specialization predicts demographic response and vulnerability of floodplain birds in Amazonia.

The annual flooding cycle of Amazonian rivers sustains the largest floodplains on Earth, which harbour a unique bird community. Recent studies suggest that habitat specialization drove different patterns of population structure and gene flow in floodplain birds. However, the lack of a direct estimate of habitat affinity prevents a proper test of its effects on population histories. In this work, we used occurrence data, satellite images and genomic data (ultra-conserved elements) from 24 bird species specialized on a variety of seasonally flooded environments to classify habitat affinities and test its influence on evolutionary histories of Amazonian floodplain birds. We demonstrate that birds with higher specialization in river islands and dynamic environments have gone through more recent demographic expansion and currently have less genetic diversity than floodplain generalist birds. Our results indicate that there is an intrinsic relationship between habitat affinity and environmental dynamics, influencing patterns of population structure, demographic history and genetic diversity. Within the floodplains, historical landscape changes have had more severe impacts on island specialists, making them more vulnerable to current and future anthropogenic changes, as those imposed by hydroelectric dams in the Amazon Basin.

Amazonian birds in more dynamic habitats have less population genetic structure and higher gene flow.

Understanding the factors that govern variation in genetic structure across species is key to the study of speciation and population genetics. Genetic structure has been linked to several aspects of life history, such as foraging strategy, habitat association, migration distance, and dispersal ability, all of which might influence dispersal and gene flow. Comparative studies of population genetic data from species with differing life histories provide opportunities to tease apart the role of dispersal in shaping gene flow and population genetic structure. Here, we examine population genetic data from sets of bird species specialized on a series of Amazonian habitat types hypothesized to filter for species with dramatically different dispersal abilities: stable upland forest, dynamic floodplain forest, and highly dynamic riverine islands. Using genome-wide markers, we show that habitat type has a significant effect on population genetic structure, with species in upland forest, floodplain forest, and riverine islands exhibiting progressively lower levels of structure. Although morphological traits used as proxies for individual-level dispersal ability did not explain this pattern, population genetic measures of gene flow are elevated in species from more dynamic riverine habitats. Our results suggest that the habitat in which a species occurs drives the degree of population genetic structuring via its impact on long-term fluctuations in levels of gene flow, with species in highly dynamic habitats having particularly elevated gene flow. These differences in genetic variation across taxa specialized in distinct habitats may lead to disparate responses to environmental change or habitat-specific diversification dynamics over evolutionary time scales.

A compreensão dos fatores que governam a variação da estrutura genética entre as espécies é fundamental para o estudo da especiação e da genética das populações. A estrutura genética tem sido ligada a vários aspectos da história da vida, tais como estratégia de forrageio, associação ao habitat, distância de migração e capacidade de dispersão, os quais poderiam influenciar a dispersão e o fluxo gênico. Estudos comparativos usando espécies que diferem nas suas histórias de vida oferecem uma oportunidade para desvendar o papel da dispersão no estabelecimento do fluxo gênico e da estrutura genética da população. Aqui examinamos dados genéticos populacionais de diversas espécies de aves com diferentes capacidades de dispersão especializadas em três habitats amazônicos, incluindo florestas de terra-firme, florestas de várzea e ilhas fluviais, cujos ambientes ripários são altamente dinâmicos. Utilizando dados genômicos que incluem milhares de loci, mostramos que o tipo de habitat tem um efeito significativo na estruturação genética das populações; espécies de florestas de terra-firme, florestas de várzea e ilhas fluviais exibem níveis de estruturação progressivamente menores. Embora os traços morfológicos frequentemente usados como indicadores da capacidade de dispersão a nível individual não expliquem este padrão, as medidas genéticas populacionais de fluxo gênico são altas em espécies associadas a habitats ribeirinhos mais dinâmicos. Nossos resultados sugerem que o habitat no qual uma espécie é encontrada determina o grau de estruturação genética da população através de seu impacto nas flutuações de longo prazo do fluxo gênico, com espécies em habitats altamente dinâmicos tendo um fluxo gênico particularmente alto. As diferenças na variação genética dos táxons especializados em diferentes habitats podem resultar em respostas díspares às mesmas mudanças ambientais, ou dinâmicas de diversificação específicas a um determinado habitat ao longo de escalas de tempo evolutivas.

Comprender los factores que rigen la variación de la estructura genética entre especies es clave para el estudio de la especiación y la genética de poblaciones. La estructura genética se ha relacionado con varios aspectos de la historia vital, como la estrategia de búsqueda de alimento, la asociación de hábitats, la distancia de migración y la capacidad de dispersión, factores todos ellos que podrían influir en la dispersión y el flujo genético. Los estudios comparativos de datos genéticos poblacionales de especies con historias vitales diferentes ofrecen la oportunidad de desentrañar el papel de la dispersión en la conformación del flujo genético y la estructura genética poblacional. En este trabajo examinamos los datos genéticos de poblaciones de especies de aves especializadas en una serie de hábitats amazónicos que, según la hipótesis, filtran especies con capacidades de dispersión radicalmente diferentes: bosques estables de tierras altas, bosques dinámicos de llanuras aluviales e islas fluviales altamente dinámicas. Utilizando marcadores genómicos, demostramos que el tipo de hábitat tiene un efecto significativo en la estructura genética de la población, y que las especies de los bosques de tierras altas, los bosques inundables y las islas fluviales presentan niveles de estructura progresivamente más bajos. Aunque los rasgos morfológicos utilizados como indicadores de la capacidad de dispersión individual no explican este patrón, las medidas genéticas poblacionales del flujo genético son más elevadas en las especies de hábitats fluviales más dinámicos. Nuestros resultados sugieren que el hábitat en el que se encuentra una especie determina el grado de estructuración genética de la población a través de su impacto en las fluctuaciones a largo plazo de los niveles de flujo genético, siendo las especies de hábitats muy dinámicos las que presentan un flujo genético particularmente elevado. Estas diferencias en la variación genética entre taxones especializados en hábitats distintos pueden dar lugar a respuestas dispares al cambio ambiental o a dinámicas de diversificación específicas del hbitat a lo largo de escalas temporales evolutivas.

Human impacts outpace natural processes in the Amazon.

Amazonian environments are being degraded by modern industrial and agricultural activities at a pace far above anything previously known, imperiling its vast biodiversity reserves and globally important ecosystem services. The most substantial threats come from regional deforestation, because of export market demands, and global climate change. The Amazon is currently perched to transition rapidly from a largely forested to a nonforested landscape. These changes are happening much too rapidly for Amazonian species, peoples, and ecosystems to respond adaptively. Policies to prevent the worst outcomes are known and must be enacted immediately. We now need political will and leadership to act on this information. To fail the Amazon is to fail the biosphere, and we fail to act at our peril.

River network rearrangements promote speciation in lowland Amazonian birds.

Large Amazonian rivers impede dispersal for many species, but lowland river networks frequently rearrange, thereby altering the location and effectiveness of river barriers through time. These rearrangements may promote biotic diversification by facilitating episodic allopatry and secondary contact among populations. We sequenced genome-wide markers to evaluate the histories of divergence and introgression in six Amazonian avian species complexes. We first tested the assumption that rivers are barriers for these taxa and found that even relatively small rivers facilitate divergence. We then tested whether species diverged with gene flow and recovered reticulate histories for all species, including one potential case of hybrid speciation. Our results support the hypothesis that river rearrangements promote speciation and reveal that many rainforest taxa are micro-endemic, unrecognized, and thus threatened with imminent extinction. We propose that Amazonian hyper-diversity originates partly from fine-scale barrier displacement processes-including river dynamics-which allow small populations to differentiate and disperse into secondary contact.

Recent population differentiation in the habitat specialist Glossy Antshrike (Aves: Thamnophilidae) across Amazonian seasonally flooded forests.

We assessed population structure and the spatio-temporal pattern of diversification in the Glossy Antshrike Sakesphorus luctuosus (Aves, Thamnophilidae) to understand the processes shaping the evolutionary history of Amazonian floodplains and address unresolved taxonomic controversies surrounding its species limits. By targeting ultraconserved elements (UCEs) from 32 specimens of S. luctuosus, we identified independent lineages and estimated their differentiation, divergence times, and migration rates. We also estimated current and past demographic histories for each recovered lineage. We found evidence confirming that S. luctuosus consists of a single species, comprising at least four populations, with some highly admixed individuals and overall similar levels of migration between populations. We confirmed the differentiation of the Araguaia River basin population (S. l. araguayae) and gathered circumstantial evidence indicating that the taxon S. hagmanni may represent a highly introgressed population between three distinct phylogroups of S. luctuosus. Divergences between populations occurred during the last 1.2 mya. Signs of population expansions were detected for populations attributed to subspecies S. l. luctuosus, but not for the S. l. araguayae population. Our results support that S. luctuosus has had a complex population history, resulting from a high dependence on southeastern "clear water" seasonally flooded habitats and their availability through time. Spatial and demographic expansions toward the western "white water" flooded forests might be related to recent changes in connectivity and availability of these habitats. Our study reinforces the view that isolation due to absence of suitable habitat has been an important driver of population differentiation within Amazonian flooded forests, but also that differences between várzeas ("white water" floodplains, mostly in southwestern Amazonia) and igapós ("clear water" floodplains, especially located in the east) should be further explored as drivers of micro-evolution for terrestrial species.

Leveraging natural history biorepositories as a global, decentralized, pathogen surveillance network.

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) pandemic reveals a major gap in global biosecurity infrastructure: a lack of publicly available biological samples representative across space, time, and taxonomic diversity. The shortfall, in this case for vertebrates, prevents accurate and rapid identification and monitoring of emerging pathogens and their reservoir host(s) and precludes extended investigation of ecological, evolutionary, and environmental associations that lead to human infection or spillover. Natural history museum biorepositories form the backbone of a critically needed, decentralized, global network for zoonotic pathogen surveillance, yet this infrastructure remains marginally developed, underutilized, underfunded, and disconnected from public health initiatives. Proactive detection and mitigation for emerging infectious diseases (EIDs) requires expanded biodiversity infrastructure and training (particularly in biodiverse and lower income countries) and new communication pipelines that connect biorepositories and biomedical communities. To this end, we highlight a novel adaptation of Project ECHO's virtual community of practice model: Museums and Emerging Pathogens in the Americas (MEPA). MEPA is a virtual network aimed at fostering communication, coordination, and collaborative problem-solving among pathogen researchers, public health officials, and biorepositories in the Americas. MEPA now acts as a model of effective international, interdisciplinary collaboration that can and should be replicated in other biodiversity hotspots. We encourage deposition of wildlife specimens and associated data with public biorepositories, regardless of original collection purpose, and urge biorepositories to embrace new specimen sources, types, and uses to maximize strategic growth and utility for EID research. Taxonomically, geographically, and temporally deep biorepository archives serve as the foundation of a proactive and increasingly predictive approach to zoonotic spillover, risk assessment, and threat mitigation.

Multiple species and deep genomic divergences despite little phenotypic differentiation in an ancient Neotropical songbird, Tunchiornis ochraceiceps (Sclater, 1860) (Aves: Vireonidae).

Several bird taxa have been recently described or elevated to full species and almost twice as many bird species than are currently recognized may exist. Defining species is one of the most basic and important issues in biological science because unknown or poorly defined species hamper subsequent studies. Here, we evaluate the species limits and evolutionary history of Tunchiornis ochraceiceps-a widespread forest songbird that occurs in the lowlands of Central America, Chocó and Amazonia-using an integrative approach that includes plumage coloration, morphometrics, vocalization and genomic data. The species has a relatively old crown age (~9 Ma) and comprises several lineages with little, if any, evidence of gene flow among them. We propose a taxonomic arrangement composed of four species, three with a plumage coloration diagnosis and one deeply divergent cryptic species. Most of the remaining lineages have variable but unfixed phenotypic characters despite their relatively old origin. This decoupling of genomic and phenotypic differentiation reveals a remarkable case of phenotypic conservatism, possibly due to strict habitat association. Lineages are geographically delimited by the main Amazonian rivers and the Andes, a pattern observed in studies of other understory upland forest Neotropical birds, although phylogenetic relationships and divergence times among populations are idiosyncratic.

Genomic phylogeography of the White-crowned Manakin Pseudopipra pipra (Aves: Pipridae) illuminates a continental-scale radiation out of the Andes.

The complex landscape history of the Neotropics has generated opportunities for population isolation and diversification that place this region among the most species-rich in the world. Detailed phylogeographic studies are required to uncover the biogeographic histories of Neotropical taxa, to identify evolutionary correlates of diversity, and to reveal patterns of genetic connectivity, disjunction, and potential differentiation among lineages from different areas of endemism. The White-crowned Manakin (Pseudopipra pipra) is a small suboscine passerine bird that is broadly distributed through the subtropical rainforests of Central America, the lower montane cloud forests of the Andes from Colombia to central Peru, the lowlands of Amazonia and the Guianas, and the Atlantic forest of southeast Brazil. Pseudopipra is currently recognized as a single, polytypic biological species. We studied the effect of the Neotropical landscape on genetic and phenotypic differentiation within this species using genomic data derived from double digest restriction site associated DNA sequencing (ddRAD), and mitochondrial DNA. Most of the genetic breakpoints we identify among populations coincide with physical barriers to gene flow previously associated with avian areas of endemism. The phylogenetic relationships among these populations imply a novel pattern of Andean origination for this group, with subsequent diversification into the Amazonian lowlands. Our analysis of genomic admixture and gene flow reveals a complex history of introgression between some western Amazonian populations. These reticulate processes confound our application of standard concatenated and coalescent phylogenetic methods and raise the question of whether a lineage in the western Napo area of endemism should be considered a hybrid species. Lastly, analysis of variation in vocal and plumage phenotypes in the context of our phylogeny supports the hypothesis that Pseudopipra is a species-complex composed of at least 8, and perhaps up to 17 distinct species which have arisen in the last ~2.5 Ma.

Phylogenomics of manakins (Aves: Pipridae) using alternative locus filtering strategies based on informativeness.

Target capture sequencing effectively generates molecular marker arrays useful for molecular systematics. These extensive data sets are advantageous where previous studies using a few loci have failed to resolve relationships confidently. Moreover, target capture is well-suited to fragmented source DNA, allowing data collection from species that lack fresh tissues. Herein we use target capture to generate data for a phylogeny of the avian family Pipridae (manakins), a group that has been the subject of many behavioral and ecological studies. Most manakin species feature lek mating systems, where males exhibit complex behavioral displays including mechanical and vocal sounds, coordinated movements of multiple males, and high speed movements. We analyzed thousands of ultraconserved element (UCE) loci along with a smaller number of coding exons and their flanking regions from all but one species of Pipridae. We examined three different methods of phylogenetic estimation (concatenation and two multispecies coalescent methods). Phylogenetic inferences using UCE data yielded strongly supported estimates of phylogeny regardless of analytical method. Exon probes had limited capability to capture sequence data and resulted in phylogeny estimates with reduced support and modest topological differences relative to the UCE trees, although these conflicts had limited support. Two genera were paraphyletic among all analyses and data sets, with Antilophia nested within Chiroxiphia and Tyranneutes nested within Neopelma. The Chiroxiphia-Antilophia clade was an exception to the generally high support we observed; the topology of this clade differed among analyses, even those based on UCE data. To further explore relationships within this group, we employed two filtering strategies to remove low-information loci. Those analyses resulted in distinct topologies, suggesting that the relationships we identified within Chiroxiphia-Antilophia should be interpreted with caution. Despite the existence of a few continuing uncertainties, our analyses resulted in a robust phylogenetic hypothesis of the family Pipridae that provides a comparative framework for future ecomorphological and behavioral studies.

Paleoclimatic evolution as the main driver of current genomic diversity in the widespread and polymorphic Neotropical songbird Arremon taciturnus.

Several factors have been proposed as drivers of species diversification in the Neotropics, including environmental heterogeneity, the development of drainage systems and historical changes in forest distribution due to climatic oscillations. Here, we investigate which drivers contributed to the evolutionary history and current patterns of diversity of a polymorphic songbird (Arremon taciturnus) that is widely distributed in Amazonian and Atlantic forests as well as in Cerrado gallery and seasonally-dry forests. We use genomic, phenotypic and habitat heterogeneity data coupled with climatic niche modelling. Results suggest the evolutionary history of the species is mainly related to paleoclimatic changes, although changes in the strength of the Amazon river as a barrier to dispersal, current habitat heterogeneity and geographic distance were also relevant. We propose an ancestral distribution in the Guyana Shield, and recent colonization of areas south of the Amazon river at ~380 to 166 kya, and expansion of the distribution to southern Amazonia, Cerrado and the Atlantic Forest. Since then, populations south of the Amazon River have been subjected to cycles of isolation and possibly secondary contact due to climatic changes that affected habitat heterogeneity and population connectivity. Most Amazonian rivers are not associated with long lasting isolation of populations, but some might act as secondary barriers, susceptible to crossing under specific climatic conditions. Morphological variation, while stable in some parts of the distribution, is not a reliable indicator of genetic structure or phylogenetic relationships.

Quaternary climate changes as speciation drivers in the Amazon floodplains.

The role of climate as a speciation driver in the Amazon has long been discussed. Phylogeographic studies have failed to recover synchronous demographic responses across taxa, although recent evidence supports the interaction between rivers and climate in promoting speciation. Most studies, however, are biased toward upland forest organisms, while other habitats are poorly explored and could hold valuable information about major historical processes. We conducted a comparative phylogenomic analysis of floodplain forest birds to explore the effects of historical environmental changes and current connectivity on population differentiation. Our findings support a similar demographic history among species complexes, indicating that the central portion of the Amazon River basin is a suture zone for taxa isolated across the main Amazonian sub-basins. Our results also suggest that changes in the fluvial landscape induced by climate variation during the Mid- and Late Pleistocene drove population isolation, leading to diversification with subsequent secondary contact.

Diversity and evolution of Amazonian birds: implications for conservation and biogeography.

Amazonia has been a focus of interest since the early days of biogeography as an intrinsically complex and extremely diverse region. This region comprises an intricate mosaic that includes diverse types of forest formations, flooded environments and open vegetation. Increased knowledge about the distribution of species in Amazonia has led to the recognition of complex biogeographic patterns. The confrontation of these biogeographic patterns with information on the geological and climatic history of the region has generated several hypotheses dedicated to explain the origin of the biological diversity. Genomic information, coupled with knowledge of Earth's history, especially the evolution of the Amazonian landscape, presents fascinating possibilities for understanding the mechanisms that govern the origin and maintenance of diversity patterns in one of the most diverse regions of the world. For this we will increasingly need more intense and coordinated interactions between researchers studying biotic diversification and the evolution of landscapes. From the interaction between these two fields of knowledge that are in full development, an increasingly detailed understanding of the historical mechanisms related to the origin of the species will surely arise.

Diversification history in the Dendrocincla fuliginosa complex (Aves: Dendrocolaptidae): Insights from broad geographic sampling.

Dendrocincla woodcreepers are ant-following birds widespread throughout tropical America. Species in the genus are widely distributed and show little phenotypic variation. Notwithstanding, several subspecies have been described, but the validity of some of these taxa and the boundaries among them have been discussed for decades. Recent genetic evidence based on limited sampling has pointed to the paraphyly of D. fuliginosa, showing that its subspecies constitute a complex that also includes D. anabatina and D. turdina. In this study we sequenced nuclear and mitochondrial markers for over two hundred individuals belonging to the D. fuliginosa complex to recover phylogenetic relationships, describe intraspecific genetic diversity and provide historical biogeographic scenarios of diversification. Our results corroborate the paraphyly of D. fuliginosa, with D. turdina and D. anabatina nested within its recognized subspecies. Recovered genetic lineages roughly match the distributions of described subspecies and congruence among phylogenetic structure, phenotypic diagnosis and distribution limits were used to discuss current systematics and taxonomy within the complex, with special attention to Northern South America. Our data suggest the origin of the complex in western Amazonia, associated with the establishment of upland forests in the area during the early Pliocene. Paleoclimatic cycles and river rearrangements during the Pleistocene could have, at different times, both facilitated dispersal across large Amazonian rivers and the Andes and isolated populations, likely playing an important role in differentiation of extant species. Previously described hybridization in the headwaters of the Tapajós river represents a secondary contact of non-sister lineages that cannot be used to test the role of the river as primary source of diversification. Based on comparisons of D. fuliginosa with closely related understory upland forest taxa, we suggest that differential habitat use could influence diversification processes in a historically changing landscape, and should be considered for proposing general mechanisms of diversification.

A dynamic continental moisture gradient drove Amazonian bird diversification.

The Amazon is the primary source of Neotropical diversity and a nexus for discussions on processes that drive biotic diversification. Biogeographers have focused on the roles of rivers and Pleistocene climate change in explaining high rates of speciation. We combine phylogeographic and niche-based paleodistributional projections for 23 upland terra firme forest bird lineages from across the Amazon to derive a new model of regional biological diversification. We found that climate-driven refugial dynamics interact with dynamic riverine barriers to produce a dominant pattern: Older lineages in the wetter western and northern parts of the Amazon gave rise to lineages in the drier southern and eastern parts. This climate/drainage basin evolution interaction links landscape dynamics with biotic diversification and explains the east-west diversity gradients across the Amazon.

Evidence for mtDNA capture in the jacamar Galbula leucogastra/chalcothorax species-complex and insights on the evolution of white-sand ecosystems in the Amazon basin.

Jacamar species occur throughout Amazonia, with most species occupying forested habitats. One species-complex, Galbula leucogastra/chalcothorax, is associated to white sand ecosystems (WSE). Previous studies of WSE bird species recovered shallow genetic structure in mtDNA coupled with signs of gene flow among WSE patches. Here, we characterize diversification of the G. leucogastra/chalcothorax species-complex with dense sampling across its distribution using mitochondrial and genomic (Ultraconserved Elements, UCEs) DNA sequences. We performed concatenated likelihood and Bayesian analysis, as well as a species-tree analysis using ∗BEAST, to establish the phylogenetic relationships among populations. The mtDNA results recovered at least six geographically-structured lineages, with G. chalcothorax embedded within lineages of G. leucogastra. In contrast, both concatenated and species-tree analyses of UCE data recovered G. chalcothorax as sister to all G. leucogastra lineages. We hypothesize that the mitochondrial genome of one of the G. leucogastra lineage (Madeira) was captured into G. chalcothorax in the past. We discuss how WSE evolution and the coevolution of mtDNA and nuclear genes might have played a role in this apparently rare event.

Connecting Amazonian, Cerrado, and Atlantic forest histories: Paraphyly, old divergences, and modern population dynamics in tyrant-manakins (Neopelma/Tyranneutes, Aves: Pipridae).

Several biogeographic hypotheses have been proposed to explain connections between Amazonian and Atlantic forest biotas. These hypotheses are related to the timing of the connections and their geographic patterns. We performed a phylogeographic investigation of Tyrant-manakins (Aves: Pipridae, Neopelma/Tyranneutes) which include species inhabiting the Amazon and Atlantic forests, as well as gallery forests of the Cerrado. Using DNA sequence data, we determined phylogenetic relationships, temporal and geographic patterns of diversification, and recent intraspecific population genetic patterns, relative to the history of these biomes. We found Neopelma to be a paraphyletic genus, as N. chrysolophum is sister to Neopelma + Tyranneutes, with an estimated divergence of approximately 18 Myrs BP, within the oldest estimated divergence times of other Amazonian and Atlantic forest avian taxa. Subsequent divergences in the group occurred from Mid Miocene to Early Pliocene and involved mainly the Amazonian species, with an expansion into and subsequent speciation in the Cerrado gallery forests by N. pallescens. We found additional structure within N. chrysocephalum and N. sulphureiventer. Analysis of recent population dynamics in N. chrysocephalum, N. sulphureiventer, and N. pallescens revealed recent demographic fluctuations and restrictions to gene flow related to environmental changes since the last glacial cycle. No genetic structure was detected across the Amazon River in N. pallescens. The tyrant-manakins represent an old historical connection between the Amazon and Atlantic Forest.

Phenotypic and Genetic Structure Support Gene Flow Generating Gene Tree Discordances in an Amazonian Floodplain Endemic Species.

Before populations become independent evolutionary lineages, the effects of micro evolutionary processes tend to generate complex scenarios of diversification that may affect phylogenetic reconstruction. Not accounting for gene flow in species tree estimates can directly impact topology, effective population sizes and branch lengths, and the resulting estimation errors are still poorly understood in wild populations. In this study, we used an integrative approach, including sequence capture of ultra-conserved elements (UCEs), mtDNA Sanger sequencing and morphological data to investigate species limits and phylogenetic relationships in face of gene flow in an Amazonian endemic species (Myrmoborus lugubris: Aves). We used commonly implemented species tree and model-based approaches to understand the potential effects of gene flow in phylogenetic reconstructions. The genetic structure observed was congruent with the four recognized subspecies of M. lugubris. Morphological and UCEs data supported the presence of a wide hybrid zone between M. l. femininus from the Madeira river and M. l. lugubris from the Middle and lower Amazon river, which were recovered as sister taxa by species tree methods. When fitting gene flow into simulated demographic models with different topologies, the best-fit model indicated these two taxa as non-sister lineages, a finding that is in agreement with the results of mitochondrial and morphological analyses. Our results demonstrated that failing to account for gene flow when estimating phylogenies at shallow divergence levels can generate topological uncertainty, which can nevertheless be statistically well supported, and that model testing approaches using simulated data can be useful tools to test alternative phylogenetic hypotheses.